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Abstract 

Laboratory  and  numerical  experiments  are  carried  out  in  a  parameter  regime 
relevant  to  open-ocean  deep  convection.  We  consider  the  case  of  convection  in  a 
rotating  stratified  ocean  of  finite  depth.  Convection  is  induced  from  the  surface  by 
extracting  buoyancy  over  a  circular  area.  The  external  parameters  are  buoyancy 
forcing  of  strength  B0  applied  over  a  circular  area  of  radius  Rs,  the  rotation  rate 
is  measured  by  /,  ambient  stratification  N  and  finite  depth  H.  Theoretical  scaling 
predictions  axe  derived  to  describe  the  length  and  velocity  scales  of  the  convective 
chimney  as  it  adjusts  under  gravity  and  rotation,  and  breaks  up  through  baroclinic 
instability.  The  scales  of  interest  include  the  number  (M),  size  (leddy)  and  strength 
(urim)  of  the  baroclinic  eddies  formed.  Also  of  interest  are  the  final  depth  of  pen¬ 
etration  of  the  convective  mixed  layer  ( d/inai )  and  the  final  volume  of  convectively 
produced  water  ( V/w )•  These  scales  are  tested  against  the  laboratory  and  numerical 
experiments  and  found  to  be  appropriate.  We  show  that  for  this  idealized  problem 
dfinai  depends  only  on  the  size  and  strength  of  the  forcing  and  the  ambient  stratifica¬ 
tion  encountered  by  the  convection  event;  it  does  not  depend  explicitly  on  rotation. 
The  volume  of  convectively  modified  water  produced  continues  to  increase  as  long 
as  forcing  continues,  but  the  rate  of  production  depends  again  only  on  the  size  and 
strength  of  the  forcing  and  on  the  ambient  stratification.  The  implications  of  the 
work  to  deep  water  formation  in  the  Labrador  Sea  and  elsewhere  are  discussed.  Fi¬ 
nally,  the  study  has  relevance  to  the  role  and  representation  of  baroclinic  eddies  in 
large-scale  circulation  of  the  ocean. 
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Chapter  1 


Introduction 


1.1  Dynamics  of  chimney  convection 

This  study  is  an  attempt  to  better  understand  the  process  of  open  ocean  deep  convec¬ 
tion.  It  is  a  complex  phenomenon  comprised  of  many  aspects  only  some  of  which  are 
well  understood.  In  this  study  we  set  up  and  study  a  highly  idealized,  but  physically 
relevant  system:  convection  in  a  stratified  rotating  fluid,  driven  by  buoyancy  loss  over 
a  circular  area  at  the  surface.  Length,  time  and  velocity  scales  for  the  key  processes 
are  deduced,  as  well  as  what  we  believe  are  the  relevant  non-dimensional  numbers 
describing  the  physics  involved.  Finally,  we  test  the  appropriateness  of  the  scalings 
by  laboratory  and  numerical  experiment. 

Open  ocean  deep  convection  is  the  phenomenon  by  which  surface  water  is 
made  dense  through  cooling  and/or  salinification  and  subsequently  sinks,  sometimes 
to  great  depth.  The  externally  imposed  parameters  are  the  buoyancy  forcing,  the 
background  stratification  of  the  fluid,  the  rotation  of  the  Earth,  and  the  total  depth  of 
the  ocean.  Buoyancy  forcing  is  characterized  by  a  magnitude  Ba  (m2/s3),  (related  for 
example  to  heat  flux  by  the  relation  Q  —  ) ,  applied  over  a  circular  area  of  radius 

R„  ( m ).  Stratification  is  measured  by  the  Brunt- Vaisala  frequency  N  =  ^ 

(s-1),  rotation  by  the  Coriolis  frequency  /  (s-1).  Finally,  ocean  depth  is  denoted 
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Figure  1.1:  A  schematic  representation  of  the  idealized  convection  model  under  con¬ 
sideration.  Buoyancy  ( Ba )  is  extracted  over  a  circular  region  at  the  surface  with 
radius  R„.  The  total  depth  is  fixed  at  H.  There  are  no  horizontal  boundaries  af¬ 
fecting  horizontal  evolution;  all  experiments  are  terminated  before  side  wall  effects 
interfere.  The  entire  system  is  under  rotation,  /. 


here  by  H  (to).  These  parameters  are  all  shown  schematically  in  Figure  1.1.  In 
addition  to  these  parameters,  it  will  be  necessary  to  consider  briefly  the  viscosity  and 
diffusivity  of  the  fluid,  v  and  k  (to2/s). 

From  these  parameters,  considered  together  with  time  t  (s),  we  can  deduce 
other  scales  for  the  features  of  a  convective  event.  The  buoyancy  forcing  drives  isolated 
overturning  cells  called  plumes,  some  of  the  details  of  which  are  described  in  the  re¬ 
view  article  [Marshall,  Whitehead  and  Yates,  1994]  and  [Jones  and  Marshall,  1993]. 
These  plumes,  perhaps  less  than  1  km  in  diameter,  act  in  concert  to  form  a  chimney 
of  well  mixed  dense  fluid  as  seen  in  Figures  1.2  through  1.5.  This  chimney  may  be 
many,  perhaps  hundreds  of  plumes  wide.  Figure  1.3  shows  a  chimney  observed  in 
the  Mediterranean  Sea.  This  event  was  observed  by  [THETIS,  1994]  during  Febru¬ 
ary  1992.  The  contours  shown  here  are  potential  density,  sigma-theta.  Figure  1.4 
shows  the  phenomenon  in  the  laboratory.  Figure  1.5  shows  a  cross  section  from  one 
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Figure  1.2:  Schematic  diagram  of  a  typical  convection  event  cross  section.  Note  the 
depth  scale  measured  from  the  surface  to  the  spreading  level. 


of  our  numerical  experiments.  Figure  1.2  shows  a  schematic  of  the  chimney  cross  sec¬ 
tion,  illustrating  our  particular  interest  in  those  processes  which  control  the  depth, 
h ,  reached  by  the  convectively  produced  mixed  layer. 

Strong  horizontal  density  gradients  are  set  up  across  the  wall  of  the  chimney 
where  isopycnals  bow  up  from  depth  to  the  sea  surface.  Under  the  influences  of  rota¬ 
tion  the  chimney  begins  to  adjust  to  geostrophic  balance.  At  the  surface,  convergence 
results  in  a  cyclonic  rim  current;  at  depth,  divergence  results  in  an  anti-cyclonic  rim 
current.  These  circulations  are  also  shown  in  Figure  1.2. 

The  opposing  currents  at  the  surface  and  at  depth  are  associated,  by  thermal 
wind,  with  horizontal  buoyancy  gradients.  The  resulting  vertical  shear  leads  to  baro- 
clinic  instability.  If  the  horizontal  extent  of  the  chimney  exceeds  the  Rossby  radius 
of  deformation  the  chimney  will  break  up  into  smaller,  Rossby-radius  scale  eddies. 
This  breakup  is  illustrated  in  Figures  1.6  through  1.9  .  Figure  1.6  shows  a  schematic 
of  the  horizontal  evolution  of  the  chimney,  illustrating  the  horizontal  scales  we  wish 
to  investigate:  uTim  ( m/s )  the  velocity  of  the  rim  current,  letuy  (m)  the  length  scale 
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Figure  1.3:  [THETIS,  1994],  cross  section  of  Meditteranean  convection  event,  seen 
using  accoustic  profiling. 
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Figure  1.4:  Photograph  of  laboratory  experiment;  shown  here  is  a  typical  cross  sec¬ 
tion.  Note  the  constriction  of  the  chimney  near  the  source,  and  the  spreading  level 
at  depth. 
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Figure  1.5:  Cross  section  plotted  from  numerical  model;  shown  here  are  the  salt 
contours.  Recall  that  temperature  is  held  constant  in  the  model,  so  that  isohalines 
are  equivalent  to  isopycnals. 

the  baroclinic  eddies,  and  M,  the  number  of  eddies  formed.  Figure  1.7  shows  a  salinity 
map  from  a  [MEDOC,  1970]  surveyed  convection  site.  Figure  1.8  shows  an  overhead 
view  of  convection  in  the  laboratory  after  the  onset  of  instability.  This  event  is  a  typ¬ 
ical  example  of  one  of  our  high  mode  laboratory  experiments.  Photographs  of  lower 
mode  instabilities  will  be  presented  below.  Figure  1.9  shows  the  same  phenomenon 
manifest  in  the  numerical  model  results. 

One  objective  of  this  thesis  is  to  determine  how  the  aforementioned  scales  de¬ 
pend  on  the  externally  imposed  parameters  and  the  non-dimensional  numbers  which 
govern  the  system.  These  non-dimensional  numbers  are  chosen  to  be  R*  —  *  ^ — , 
the  natural  Rossby  number  which  measures  the  significance  of  the  distance  over  which 
rotation  becomes  important,  versus  the  total  depth  of  the  fluid  -  for  a  review  see, 
[Marshall,  Whitehead  and  Yates,  1994],  or  [Maxworthy  and  Narimousa,  1994],  N/f, 
relating  stratification  and  rotation,  and  Rs/H,  an  aspect  ratio  comparing  the  size  of 
the  forcing  area  to  the  total  depth  of  the  fluid. 
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Figure  1.6:  Schematic  diagram  of  a  typical  convection  event  as  seen  from  above, 
illustrating  our  interest  in  the  number  of  eddies  formed,  their  associated  velocities 
and  length  scales. 


Figure  1.7:  Salinity  map  showing  a  convection  region  in  the  Meditteranean  Sea, 
surveyed  hydrographically  by  [MEDOC,  1970].  Note  the  geometry  and  the  horizontal 
scale  of  the  convective  patch  -  approximately  100  km. 
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Figure  1.8:  Photograph  of  laboratory  experiment;  shown  here  is  a  typical  high  mode 
laboratory  experiment.  Note  particularly  the  distinct  eddies  migrating  away  from  the 
source. 


Figure  1.9:  Horizontal  map  of  numerical  model  velocity  fields.  Note  the  developed 
eddy  circulations,  allowing  estimates  of  mode  number,  size  and  associated  velocity. 
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To  test  the  proposed  scalings  of  h,  uTim,  leddy  and  M,  and  determine  the  de¬ 
pendence  of  these  scales  on  the  non-dimensional  numbers,  we  carried  out  laboratory 
and  numerical  experiments  simulating  our  idealized  physical  system.  The  laboratory 
experiment  consists  of  a  rotating  tank  filled  with  linearly  stratified  fluid  from  which 
buoyancy  is  extracted  over  a  circular  area  at  the  surface.  The  computational  analog 
is  a  numerical  simulation  of  the  laboratory  experiment  using  a  non-hydrostatic  prim¬ 
itive  equation  model  [Brugge  et  al.,  1991].  The  study  also  provides  an  interesting 
context  in  which  to  compare  in  detail  laboratory  and  numerical  results. 

In  Chapter  2  we  introduce  non-dimensional  numbers  and  the  proposed  scalings 
for  h ,  urim,  leddy  and  M.  Chapter  3  presents  the  laboratory  experiment,  including 
method  and  results.  Chapter  4  presents  the  numerical  simulations  of  the  laboratory 
experiments,  including  method,  comparison  with  the  laboratory,  and  further  results. 
Chapter  5  is  a  discussion  of  our  results  in  the  context  of  the  ocean. 
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Chapter  2 


Scaling 


2.1  Non-dimensional  numbers 


When  an  unstratified  fluid  of  total  depth  H  is  forced  to  convect  by  uniform  surface 
cooling  in  the  presence  of  rotation  (/),  there  is  only  one  nondimensional  combination 
of  the  main  external  parameters  Ba,  f  and  H:  the  natural  Rossby  number  R*a. 


K  ~ 


sfBM*  Ire 

— -  - 

H  H' 


Here  lTOt  =  y  B0j  f3  is  the  vertical  scale  to  which  convection  penetrates  in 
an  unstratified  fluid  during  an  inertial  period  -  the  period  when  rotation  becomes 
important.  If  R *  is  large  then  the  convection  will  be  limited  by  the  ocean  depth 
before  it  feels  the  effects  of  rotation;  if  R*a  is  small  then  the  convection  will  come 
under  geostrophic  control  before  encountering  the  bottom.  Oceanic  values  of  R*  in 
convective  regions  are  typically  small,  placing  the  phenomenon  in  a  regime  where 
rotation  may  be  important.  In  the  ocean,  typical  values  for  f,  B0  and  H  are:  in  the 
Mediterranean  /  «  10-45_1,  Ba  ~  10 ~7m2/s3  (corresponding  to  a  heat  loss  of  ~  2000 
W/m2,  where  Ba  =  a  js  the  thermal  expansion  coefficient  for  water,  p  is  the 
density  of  water,  Cp  is  the  specific  heat,  g  is  the  gravitational  constant  and  Q  is  the 
heat  loss,  and  H  ~  2000m  [Leaman  and  Schott]  giving  an  R*  of  0.1;  in  the  polax 
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seas  these  values  are  /  sa  1.5  x  10_4s_1,  Ba  ~  10 ~7m2/s3,  and  H  ~  1000  —  4000m 
[Marshall,  Whitehead  and  Yates,  1994]  giving  an  R*a  of  0.01.  Typical  ocean  values 
then  are  in  the  range  of  0.01  <  R*  <  0.1.  In  the  laboratory  experiments  presented 
here  R*a  ranges  from  0.04-0.41.  In  our  numerical  experiments  R*  varies  from  0.05-0.58. 
In  both  cases  the  experiments  span  the  oceanic  regime. 

In  addition  to  R*a  two  additional  non-dimensional  numbers  are  needed  to  fully 
describe  the  system  when  stratification  and  a  finite  forcing  length  scale  are  included, 

N/f 

and 

Rs/H. 

The  first  describes  the  relative  significance  of  rotation  and  stratification.  The 
second  is  an  aspect  ratio  describing  the  size  of  the  forcing  area  compared  to  the 
depth  of  the  water  column.  Values  of  N  for  the  upper  ocean  are  typically  0.01  s_1 
in  the  high  gradient  regions,  while  values  for  the  deep  ocean  may  be  as  small  as 
0.001  s_1,  [Gill,  1982].  Typical  values  of  /  are  of  order  10-4s-1.  Thus  the  range 
of  N/f  for  the  ocean,  specifically  for  convection  regions  is  from  order  1-10,  while 
values  for  the  rest  of  the  ocean  may  be  as  high  as  N/f  ~  50.  Deep  convection 
sites  are  very  weakly  stratified,  [Leaman  and  Schott].  In  the  laboratory  we  achieve 
a  range  of  0.86  <  N/f  <  6.64;  in  the  numerical  experiments  we  achieve  a  range  of 
0.82  <  N/f  <  7.5.  This  places  the  experiments  described  here  in  the  lower  limit  of 
thermocline  ocean  values,  making  the  experiments  particularly  suited  to  deep  ocean 
convection  sites.  Note  also  that  N/f  spans  an  order  of  magnitude  in  the  experiments 
to  allow  investigation  of  dependence  of  convection  on  this  parameter.  Typical  length 
scales  for  R,  are  set  by  meteorological  forcing  and  in  the  ocean  can  range  from  small 
disturbance  scales  (eg.  squall  line)  .5  km,  to  atmospheric  meso-scale  disturbances 
100  km.  Ocean  values  of  H  are  1000-4000  m.  This  gives  a  range  of  .1  <  Ra/H  <  100. 
In  both  the  laboratory  and  the  numerical  experiments  Ra/H  spanned  .25-. 75.  This 
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places  the  experiments  again  at  the  lower  limit  of  the  real  ocean  values,  but  still  well 
within  the  range  that  makes  them  relevant  to  the  real  ocean. 

One  other  non-dimensional  number  must  be  mentioned,  that  is  the  Rayleigh 
number,  Rdf  =  It  relates  the  relevant  forcing  and  stabilizing  mechanisms 

involved  in  convective  overturning.  Buoyancy  forcing  at  the  surface  establishes  a 
density  gradient;  diffusion  erodes  that  gradient.  It  is  a  well  established  result  from 
Rayleigh’s  calculations  that  for  a  given  stratification,  diffusion,  viscosity  and  depth 
there  exists  a  critical  value  of  forcing  necessary  to  achieve  overturning.  In  other 
words,  there  is  a  critical  value  of  Rdf  necessary  to  achieve  instability.  For  the  simple 
problem  of  temperature  difference  between  two  plates,  that  number  is 

Rdf(cr)  =  — ~ —  =  656. 

This  result  is  highly  dependent  on  the  boundary  conditions.  The  result  will  also  be  af¬ 
fected  by  the  presence  of  rotation.  Rotation  stiffens  the  fluid  and  inhibits  overturning, 
then  the  critical  value  of  Rdf  will  be  increased. 

The  critical  value  has  been  related  to  the  Taylor  number  [Chandrasekhar,  1953] 
and  [Nakagawa  and  Frenzen,  1955]  according  to  the  relationship 

Rdf(cr )  =  kTd 2/3, 

where  k  is  a  constant  dependent  on  the  boundary  conditions,  and  Td  is  the  Taylor 
number  Td  =  . 

[Bubnov  and  Senatorsky,  1988]  determined  k  experimentally  and  found  it  to 
lie  between  2.39  and  8.72;  the  particular  value  depended  strongly  on  the  boundary 
conditions  used.  For  the  laboratory  and  numerical  experiments  studied  here,  Td  ^ 
105  giving  a  critical  Rayleigh  number  of  Rdf(cr)  =  106. 

The  important  point  to  be  made  is  that  there  is  a  critical  value  of  Rdf  necessary 
for  convective  overturning  and  that  number  Rdf(cr )  is  never  larger  than  106  at  the 
upper  limit.  In  our  numerical  and  laboratory  experiments  we  exceed  that  number  by 
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several  orders  of  magnitude.  In  the  ocean,  /  ta  10  4/s,  H  m  103m,  v  «  10  6m2/s 
and  k  ~  10~7m2/s  (for  thermal  diffusion  in  water).  This  gives  typical  ocean  values 
of  Rdf  «  1027  and  To  «  1022. 

In  the  laboratory  it  is  impossible  to  achieve  a  value  of  H  similar  to  that  of  the 
ocean.  Ocean  values  of  H  may  be  obtained  in  the  numerical  experiments,  but  only  at 
the  price  of  increasing  values  of  v  and  k.  The  increase  in  H  is  accompanied  by  a  loss 
of  resolution  which  must  be  compensated  by  increases  in  v  and  k  to  parameterize  all 
unresolved  mixing.  So,  using  molecular  values  of  v  and  k,  typical  laboratory  values 
of  Rdf  range  from  1013  to  1015  where  values  of  /  range  from  .01  to  .1/s  and  H  ~  ,1m 
resulting  in  Rdf. 

In  the  numerical  experiments  a  larger  values  of  k.  is  used  for  numerical  stability, 
k  fa  10 ~6m2/s.  So,  the  numerical  system  has  a  Prandtl  number  of  one.  Rotation 
can  be  increased  more  in  the  numerical  experiments  than  in  the  laboratory,  so  /  fa 
.01  —  1  /s.  Values  of  Rdf  for  the  numerical  experiments  range  from  107  —  1011. 

In  summary,  it  is  impossible  to  achieve  oceanic  values  of  Rdf  in  either  the 
laboratory  or  the  numerical  experiments.  However,  it  is  still  possible  to  gain  useful 
information  from  the  experiments  because  in  both  the  laboratory  and  on  the  com¬ 
puter,  the  critical  value  of  Rdf  is  exceeded  by  several  orders  of  magnitude,  placing 
all  three  cases  in  the  realm  of  convective  instability.  So,  we  can  use  laboratory  and 
numerical  experiments  to  systematically  investigate  convective  instability  and  gain 
meaningful  results  which  may  be  applied  to  ocean  dynamics. 


2.2  Convection  in  a  stratified  rotating  ocean 

We  have  chosen  to  study  the  problem  of  convection  in  a  stratified  rotating  ocean  with 
a  finite  depth.  We  will  consider  buoyancy  forcing  applied  over  a  finite  circular  area 
only.  In  our  highly  idealized  problem  the  parameters  to  be  considered  are:  B0,  t,  N, 
f,  H,  and  R„.  Here  again  Ba  ( m2/s 3)  is  the  magnitude  of  the  buoyancy  forcing,  t  (s) 
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is  time,  N  (s-1)  is  the  Brunt- Vaisala  frequency  and  /  (s-1)  the  Coriolis  frequency, 
H(m )  is  the  total  depth  of  the  ocean,  and  Rs  is  the  radius  of  the  buoyancy  forcing 
patch. 

Previous  work  [Legg  and  Marshall,  1994]  [Jones  and  Marshall,  1993]  has  fo¬ 
cused  on  the  evolution  of  individual  plumes.  In  the  present  study  we  concern  ourselves 
with  the  development  of  the  entire  chimney  which  is  made  up  of,  and  indeed  created 
by,  many  plumes.  The  plumes  act  in  concert  to  densify  the  surface  water  layer.  The 
entire  area  of  the  chimney  is  mixed,  the  dense  layer  evolving  with  time.  One  premise 
of  this  scaling  is  that  the  deepening  of  the  convective  mixed  layer  is  independent  of  the 
inhomogeneity  introduced  by  the  existence  of  the  individual  plumes.  The  assumption 
that  the  chimney  evolves  independent  of  the  resolution  of  the  plumes  is  supported  by 
both  laboratory  and  computer  data  and  is  discussed  in  Chapter  5.  The  chimney,  if  it 
is  sufficiently  large,  will  ultimately  break  up  under  baroclinic  instability.  Recall  that 
the  properties  of  the  convective  chimney  which  we  have  chosen  to  look  at  are  h,  uTim, 
leddy  and  M. 

Let  us  now  consider  the  vertical  evolution  of  the  mixed  layer  at  a  position  in 
the  center  of  the  well  mixed  chimney.  A  water  parcel  in  this  position  will  be  unaware 
of  any  spatial  inhomogeneity  in  either  cooling  or  overturning. 

Before  rotation  becomes  important  the  evolution  of  the  mixed  layer  can  be 
described  by  a  one  dimensional  model,  assuming  there  is  no  lateral  mixing  across  the 
boundary  of  the  chimney.  From  the  definition  of  heat  flux,  one  can  write 


(  dT\ 


where  h  is  the  depth  over  which  the  fluid  is  mixed,  T  is  temperature,  and  t  is  time  and 
the  time  derivative  of  temperature  has  been  written  thus:  =  fjff-  The  resulting 

expression  for  buoyancy  is: 


B0  =  —goth 


17 


Figure  2.1:  A  schematic  diagram  of  the  one  dimensional  vertical  evolution  of  a  mixed 
layer  driven  by  temperature  forcing.  Each  graph  represents  a  new  time.  As  buoyancy 
is  removed  through  the  surface,  the  mixed  layer  extends  to  a  larger  depth. 


where  Ba  — 

pep 


Grouping  terms  and  using  the  definition  of  buoyancy  frequency  N, 


AT2  9T 

N  =  ~9a~dz ’ 


we  find 


B0  =  N2h 


dh 


Integrating  and  solving  this  equation  for  h  gives  the  well-known  result  for  the  non- 
penetrative  deepening  of  a  mixed  layer  (see  for  example  [Ivey  et  al.,  1993]): 


h  = 


This  derivation  assumes  only  that  the  mixing  is  non-penetrative.  The  result 
suggests,  as  expected,  that  if  N  increases  then  h  decreases;  that  is,  larger  stratification 
will  impede  the  deepening  of  the  mixed  layer.  It  also  suggests  that  for  very  small 
times  h  will  increase  rapidly,  while  for  large  times  h  increases  much  more  slowly. 
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We  might  expect  from  geostrophic  adjustment  theory  that  the  chimney  wall 
has  a  width  Ru  =  the  Rossby  radius  of  deformation.  As  rotation  becomes 
important  the  chimney  adjusts  toward  geostrophic  balance.  The  surface  convergence 
described  in  Chapter  1  has  an  associated  cyclonic  velocity,  Using  thermal  wind 

balance  one  can  deduce  uTim: 

du  g  ( dp\ 

dz  pf  V&/ j  ’ 

where  u  is  the  azimuthal  velocity  and  indicates  the  density  gradient  across  the 
chimney  wall.  Then 

gh  Apy 

UTim  —  r  a 

Pol  Ay 

We  define  the  chimney  geometry  such  that  the  density  difference  radially  out¬ 
ward  from  the  center  is  the  same  as  the  vertical  density  difference  from  the  spreading 
level  to  the  surface,  see  Figure  2.2.  Then  by  definition  Apy  =  A pz.  So 


Substituting  into  our  expression  for  ur;m 

hgApzhL  N2h2 
Urim  ~  fPoAyhL  ~  fL 

which,  recalling  that  L  =  gives 


~Nh. 


Thus  the  rim  current  will  continue  to  grow  as  h  evolves  in  time,  until  the  onset 
of  baroclinic  instability.  At  that  time  the  eddies  formed  will  transport  water  away 
from  the  chimney,  violating  the  1-d  assumption  of  no  lateral  mixing. 

The  breakup  time  scale  may  be  derived  from  the  Eady  result  for  baroclinic 
instability  growth  rate.  The  growth  rate  of  the  fastest  growing  mode  is  described  by 
the  complex  part  of  the  wave  frequency;  the  result  [Eady,  1949]  being 

,  / 
kCi~VK 
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Figure  2.2:  A  schematic  cross  section  through  a  chimney,  illustrating  the  definitions 
of  chimney  geometry.  Note  especially  the  equivalence  of  the  density  differences  from 
the  center  of  the  chimney  to  the  outside,  and  from  the  spreading  layer  to  the  surface. 


where  R,  =  N2/  (g)’. 

Once  again  assuming  thermal  wind  balance  we  can  use  the  above  result  of 


g  ~  N,  ie.  the  Ri  =  1,  then 


kci  =  /, 


leading  to  a  time  scale  t  ~  /  1 . 

At  this  time,  t  ~  f~l  our  derived  scales  become: 


so 

Urim  ~  Nh  ~  y/2B0/f, 

and 

L-y/WJP. 

Interestingly,  each  of  these  quantities  appears  in  terms  of  the  fundamental 
scales  governing  the  evolution  of  the  individual  plumes  in  an  unstratified  rotating 
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fluid  [Jones  and  Marshall,  1993]: 


Irot  ~  y/B0/P 

'U'TOt  ^  \J B0j  f . 


Then 


h  ~ 


l^rot 

If’ 


l^rim  ~  'U'roti 


and 


L  ~  lTOt- 


2.3  The  equilibrium  state 

So  the  initial  stages  of  the  vertical  evolution  of  the  mixed  layer  are  governed  by 
the  one  dimensional  result,  which  is  valid  until  baroclinic  instability  sets  in  on  a 
time-scale  t  ~  /-1 .  At  that  time,  the  baroclinic  instability  drives  exchange  between 
the  chimney  and  the  surrounding  region.  This  lateral  exchange  inhibits  the  further 
vertical  evolution  of  the  mixed  layer,  so  that  the  depth  with  time  curve  now  lies 
below  the  t1?2  curve,  see  Figure  2.3.  Eventually  an  equilibrium  must  be  reached 
between  the  density  forcing  through  the  surface  and  the  lateral  removal  of  dense  fluid 
by  the  baroclinic  eddies,  as  detailed  in  [Legg  and  Marshall,  1994].  This  equilibrium 
necessarily  halts  further  deepening  of  the  mixed  layer;  the  depth  at  this  time  will 
be  the  maximum  depth  achieved  and  may  derived  as  follows.  The  balance  between 
the  buoyancy  flux  through  the  surface  and  the  lateral  transport  can  be  written  as 
[Visbeck  and  Marshall,  1994] 

/  —dA  =  -  f  i^dldZ, 

J  g  p  J  J 

where  u  is  the  horizontal  cross  stream  velocity  and  the  bar  denotes  a  time  average, 
indicating  our  desire  to  study  the  steady  state.  From  energy  analysis  of  the  baroclinic 
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region  we  can  write  [Eady,  1949]  [Green,  1970] 


"7  "T  9  12 

uy  =  WS  ' 


where  p'  is  the  perturbation  density  anomaly  achieved  when  moving  a  water  parcel 
in  the  baroclinic  region.  Since  (as  discussed  above) 


u'p'  may  be  written  thus: 


=  N2h, 


u'p'  =  —N3h. 


Then  our  above  integral  balance  becomes 


which  reduces  to 


-B0ttR2s 


h  ~  d  final  ~ 


-2itRsN3h3, 


( BoR,)1 /3 


This  is  the  quantity  we  call  dfinai,  our  scaling  for  the  maximum  depth  achieved  by  the 
convective  mixed  layer,  assuming  all  available  potential  energy  is  converted  to  kinetic 
energy.  The  resulting  scaling  for  penetration  depth  says  that  d/i„ai  should  depend 
only  on  the  size  of  the  source,  the  strength  of  the  forcing  and  stratification.  There 
is  no  /  dependence.  This  prediction  is  strongly  supported  by  the  numerical  data  in 
Chapter  4. 

Now  let  us  go  back  and  investigate  the  scales  of  the  baroclinic  instability,  in 
the  context  of  the  evolution  of  the  mixed  layer.  We  know  from  the  solution  of  Eady’s 


baroclinic  instability  problem  that 


_  Nh 
Rd-j-. 


As  h  increases  then  so  must  leddy,  in.  accordance  with  the  expressions  derived  above  for 
h.  The  number  of  eddies  present  is  dependent  on  the  integral  number  of  wavelengths 
that  may  be  contained  about  the  circumference  of  the  chimney,  that  is 

. .  2t tR2 
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^  on  ^  off  ^  instability 

Figure  2.3:  Plot  illustrating  the  arrest  of  the  vertical  evolution  of  the  mixed  layer  in 
the  presence  of  baroclinic  eddy  transport.  As  the  eddies  move  dense  water  away  from 
the  source  area,  the  mixed  layer  ceases  to  deepen.  This  is  manifest  in  a  leveling  off  of 
the  average  salt  content  in  a  cylinder  defined  by  the  source  area  and  the  total  depth 
of  the  water  column. 


Again,  this  number  will  change  with  time,  as  Ze<My  grows.  This  prediction  is  borne 
out  by  the  data  as  initially  h  is  small  and  the  chimney  begins  to  break  up  into  a 
large  number  of  eddies.  As  time  goes  on,  h  increases  and  one  sees  the  eddies  coalesce. 
In  the  end,  the  number  of  eddies  actually  shed  is  smaller  than  the  number  of  waves 
initially  visible  at  the  onset  of  instability.  In  fact,  according  to  the  analysis  of  the 
problem  for  a  cylindrical  geometry,  carried  out  in  [Legg  and  Marshall,  1994],  the  final 
result  should  be 

M  = 

where  h  is  chosen  appropriately. 

To  summarize  we  now  have  scalings  for  the  velocity  of  the  rim  current  in  the 
baroclinic  zone,  the  number  of  baroclinic  eddies  formed,  their  associated  length  scale 
and  the  final  resting  depth  of  the  convectively  produced  dense  water.  These  scaling 
predictions  are  hence  forth  denoted  by  the  superscript  ‘seal1: 


VsR2, 

Nh/f 


seal 


=  Nh 
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and 


iscal  _ 

eddy 


Nh 

T 


Mscal 


V*R2J 

Nh 


d 


seal 

final 


(RsBq)1'3 

N 


Finally,  it  is  interesting  to  consider  the  final  volume  of  dense  water  produced 
during  a  convection  event.  An  upper  limit  on  this  number  is  given  by  the  fastest  rate 
at  which  dense  water  may  be  produced.  The  shortest  time  in  which  the  system  can 
reach  dfinai  is  given  by  the  time  it  would  take  the  one-dimensional  model  to  evolve 
to  dfinai.  Once  the  one- dimensional  model  ceases  to  be  valid  vertical  evolution  has 
slowed  under  the  stiffness  imposed  by  rotation.  So,  the  shortest  time  necessary  for 
the  system  to  reach  dfinai  is  given  by 


h  ~  dfinai 

y/2 BJ  ( B0RS f>3 

-  rsj  - 

N  N 


Solving  this  for  t  gives  us: 


Vt  ~  b;1/6r)/3 


Then  the  fastest  rate  at  which  dense  water  can  be  produced  is  the  volume  of  the 
cylinder  beneath  the  source,  divided  by  this  time: 


v  at  6Tn(1T 


7 T  Rs  dfinai 

TfT' 


rater, 


7T  R^dfinnlBV3 


Then  the  largest  amount  of  water  produced  by  the  event  would  be  the  fastest  rate 
multiplied  by  the  amount  of  time  gone  by;  obviously  the  amount  of  water  convectively 
produced  will  continue  to  increase  as  long  as  the  forcing  continues.  So,  we  can  imagine 
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that  the  cylinder  will  fill,  no  faster  than  the  rate  we  calculated,  then  the  dense  water 
will  be  transported  away  and  the  process  will  begin  again.  We  can  expect  the  final 
volume  of  convectively  produced  water  to  follow  the  relation: 


T/»caI 
' final 


TQit&max  ♦  t . 


This  hypothesis  is  supported  by  the  numerical  data. 
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Chapter  3 


The  Laboratory  Analog 


3.1  Instrumental  method 

The  laboratory  apparatus  was  designed  by  Jack  Whitehead,  Bob  Frazel  and  the  au¬ 
thor.  The  raw  data  were  extracted  from  the  tapes  by  Jack  Whitehead.  The  statistical 
analysis  was  performed  by  the  author.  A  square  glass  tank  was  used  of  dimension 
114  cm  x  114  cm  x  60  cm,  see  Figure  3.1.  The  bottom  of  the  tank  is  marine  plywood. 
Dense  fluid  is  introduced  at  the  surface  through  a  diffuser.  Two  types  of  diffuser  are 
used  depending  on  the  source  area  required  for  a  given  run.  Two  video  cameras  are 
used  to  record  the  experiments,  one  at  the  side  of  the  tank,  the  other  rotating  with 
the  table,  recording  from  above. 

The  tank  is  filled  using  a  two  reservoir  method  to  achieve  stratification.  A 
fresh  water  reservoir  is  filled  with  one  half  the  volume  of  the  experimental  tank.  A 
second  reservoir  is  also  filled  with  one  half  the  tank  volume  but  has  a  density  equal 
to  the  maximum  density  desired  for  the  experimental  tank.  A  valve  connects  the  two 
reservoirs  at  their  base.  Water  is  drawn  from  the  fresh  reservoir  and  replaced  from 
the  salty  reservoir.  The  result  is  linearly  increasingly  dense  layers  being  set  down  in 
the  tank,  one  beneath  the  other.  The  experimental  tank  is  filled  while  it  is  turning, 


26 


1  camera 


tank 


source  diffuser 


VCR's  and  monitors 


*||g 

Sl?tl 


source  water 


rotating  table 


slip  rings 


Figure  3.1:  A  side  view  schematic  of  the  laboratory  apparatus  used  for  this  study. 

through  a  diffuser  on  the  bottom  covered  with  an  inch  of  gravel.  This  helps  maintain 
solid-body  rotation.  The  tank  is  filled  to  30  cm  above  the  gravel. 

The  dense  water  introduced  at  the  surface  is  marked  with  neutrally  buoyant 
blue  dye  to  facilitate  flow  visualization.  The  first  source  used  was  a  paint  sprayer; 
the  configuration  is  depicted  in  Figure  3.2.  The  sprayer  is  housed  in  a  plastic  cage 
which  serves  several  purposes.  The  cage  permits  mist  to  fall  only  on  the  desired  patch 
and  catches  any  drips  which  accumulate  on  the  sides  (ensuring  spatial  uniformity  of 
coverage).  The  cage  also  halts  any  wind  generated  by  the  paint  sprayer,  preventing 
the  transfer  of  momentum  to  or  evaporative  cooling  from  the  surface.  The  paint 
sprayer  generates  a  fine-enough  mist  to  satisfy  the  uniformity  and  zero  momentum 
criterion. 

The  second  source  used  allows  a  larger  area  to  be  affected.  This  source  is  a 
variation  on  that  used  by  [Maxworthy  and  Narimousa,  1994],  A  hard  plastic  cylinder 
of  the  desired  diameter  is  used.  The  bottom  is  covered  with  a  layer  of  mesh  topped 
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Figure  3.2:  A  schematic  of  the  paint  sprayer  used  to  introduce  dense  water  to  the 
experiments  discussed  here. 

by  a  layer  of  filter  paper,  topped  by  a  layer  of  sponge.  Dense  fluid  is  introduced  to 
this  diffuser  via  a  small  tube  leading  from  a  reservoir. 

The  laboratory  procedure  begins  with  the  tank  rotating.  It  is  filled  from  the 
bottom  with  a  stratified  fluid  and  left  undisturbed  until  it  has  reached  solid  body 
rotation  (this  is  approximately  2  hours).  When  this  is  established  the  experiment 
begins  with  the  initiation  of  the  buoyancy  flux.  The  resulting  flow  is  recorded  using 
the  video  cameras.  Surface  currents  are  marked  with  passive  drifters.  The  experiment 
terminates  when  the  flow  reaches  a  side  wall.  In  some  instances  it  was  necessary  to 
terminate  the  laboratory  experiment  sooner  due  to  a  practical  problem  with  the  paint 
sprayer  source.  To  prevent  salt  build-up  in  the  electronics  of  the  laboratory,  the  paint 
spray  source  is  always  discontinued  after  60-70  s. 

Salt  is  used,  rather  than  temperature,  to  induce  stratification  and  buoyancy 
loss.  This  eliminates  many  practical  laboratory  problems  including:  secondary  con¬ 
vection  during  filling  and  run;  unwanted  circulation  due  to  the  temperature  gradient 
between  the  center  of  the  tank  and  the  side  walls;  and  unmeasurable  buoyancy  loss 
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due  to  uncontrolled  evaporation.  Salt  diffuses  one  hundred  times  more  slowly  than 
heat  and  hence  a  salt  anomaly  is  easier  to  manipulate  than  a  temperature  anomaly. 
However,  it  was  not  possible  to  make  the  paint  spray  cage  completely  air  tight  so  salt 
mist  eventually  escapes  into  the  laboratory. 


3.2  Measurement  and  qualitative  results 

The  laboratory  experiments  were  set  up  primarily  to  explore  the  widest  possible 
range  of  N ,  /,  and  R„.  This  distribution  of  experimental  parameters  facilitates  the 
investigation  of  the  scalings  and  df^al  presented  in  Chapter  2.  Table  3.1  lists  the 
experimental  parameters  used  in  each  run.  Also  displayed  are  the  non-dimensional 
numbers  N/f,  R*0,  and  Ra/H.  In  all  cases  H  =  30  cm  was  the  total  depth  of  the  fluid 
in  the  laboratory  tank.  The  first  set  of  runs  were  performed  with  a  constant  value  of 
Rs  =  15  cm,  values  of  N  ranging  from  0.2  <  N  <  0.9  /s  and  values  of  /  ranging 
from  0.125  <  /  <  0.5  /s.  The  second  set  of  experiments  were  performed  varying 
the  value  of  R„  between  7.5  <  R,  <  22.5  cm,  again  varying  N  and  /. 

The  properties  measured  in  the  laboratory  are  the  mode  number  of  the  baro- 
clinic  instability  ( Mlab ),  length  scale  of  the  baroclinic  eddies  (^L),  the  associated 
maximum  velocity  (*4im) >  and  the  final  depth  of  convective  penetration,  {d!^bax).  The 
mode  number  is  measured  in  the  laboratory  by  observing  the  shape  of  the  chimney 
from  above  once  eddies  are  fully  formed.  Figure  3.3  (a)  shows  the  addition  of  source 
water  through  the  plexi-glass  diffuser,  as  seen  from  above.  Immediately  one  sees  one 
drawback  in  the  laboratory  set-up:  no  data  can  be  taken  from  the  area  directly  be¬ 
neath  the  source.  It  is  impossible  to  trace  velocities  when  the  source  is  obscuring  the 
view.  In  the  case  of  the  plexi-glass  source,  this  area  is  obscured  for  the  duration  of 
the  experiment.  In  the  case  of  the  paint  spray  source,  the  entire  surface  of  the  tank 
is  obscured,  but  only  up  until  the  source  is  switched  off  at  60-70  s.  Thereafter,  the 
cage  is  removed  offering  an  unimpeded  view  of  the  surface  flow. 
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Lab  run 

/(!/*) 

N(l/s) 

R,  (cm) 

N/f 

K 

Rsj H  ( H  =  30cm) 

1 

.82 

15 

0.50 

2 

Em 

.45 

15 

0.90 

.05 

0.50 

3 

15 

0.40 

.05 

0.50 

4 

0.250 

15 

3.31 

.15 

0.50 

5 

0.250 

.41 

15 

1.65 

.15 

0.50 

6 

jfgspH 

.22 

15 

0.88 

.15 

0.50 

7 

B1S5M 

15 

6.64 

.41 

0.50 

8 

15 

3.28 

m 

9 

.24 

15 

1.92 

■a 

.8 

7.5 

1.60 

j£j£] 

0.25 

11 

7.5 

.05 

0.25 

12 

.29 

7.5 

0.25 

13 

.86 

7.5 

0.25 

14 

.4 

1.60 

0.25 

15 

1MB 

.22 

7.5 

0.88 

0.25 

16 

.36 

15 

0.72 

.04 

0.50 

17 

.4 

15 

1.60 

.11 

0.50 

18 

.37 

15 

1.48 

.31 

0.50 

19 

15 

0.50 

20 

.43 

21 

■0 

0.75 

0.250 

.44 

0.75 

0.125 

.44 

.41 

0.75 

Table  3.1:  This  table  displays  the  experimental  parameters  for  the  laboratory  exper¬ 
iments  reported  here.  Also  shown  are  the  values  of  the  significant  non-dimensional 
parameters.  The  first  9  experiments  listed  are  those  used  in  comparison  with  the 
numerical  estimates.  The  rest  are  used  to  further  test  the  scaling  theory  for  dfC£al. 
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Figure  3.3:  Photographs  of  the  laboratory  experiment:  a)  overhead  view  ot  laboratory 
experiment  at  start,  note  the  dyed  water  just  becoming  apparent  at  edges  of  source, 
and  reservoir  of  dyed  fluid  present  in  the  source  diffuser;  b)  side  view  of  the  laboratory 
experiment,  note  the  narrow  chimney  beneath  the  source,  and  the  spreading  of  the 
dyed  fluid  at  depth;  c)  overhead  view  of  laboratory  experiment,  two  distinct  eddies 
have  formed  and  are  propagating  away  from  the  source;  d)  again  from  overhead,  now 
the  anti-cyclonic  circulation  is  accreting  dyed  water  from  the  cyclonic  eddy. 
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Once  the  dense  chimney  fluid  reaches  its  neutrally  buoyant  depth,  it  begins  to 
spread  out  from  under  the  source  and  becomes  visible  from  the  above.  The  horizontal 
extent  of  the  spreading  is  evident  in  the  cross  section  photograph  displayed  in  Figure 
3.3  (b).  The  chimney  breaks  into  a  discrete  number  of  eddies  which  migrate  away 
from  the  source.  In  many  instances,  new  eddies  seen  from  above  form  under  the 
source  when  the  first  eddies  drift  away.  The  number  of  eddies  being  simultaneously 
shed  by  the  chimney  is  identified  as  the  mode  number.  As  discussed  above,  the  initial 
mode  number  that  appears  in  the  experiment  decreases,  as  the  chimney  penetrates 
deeper,  the  visible  effect  being  a  coalescence  of  eddies.  Figure  3.3  (c)  shows  two  well 
developed  eddies  propagating  away  from  the  source.  They  are  each  still  attached  to 
the  source  water  through  a  narrow  filament.  The  following  panel  shows  the  two  eddies 
at  a  later  time  when  some  of  the  dye  is  being  advected  around  into  the  anti-cyclone. 

Recall  that  Figure  1.8  showed  a  high  mode  instability,  illustrating  a  surface 
manifestation  of  the  sub-surface  anti-cyclone.  The  surface  cyclone  and  sub-surface 
anti-cyclone  always  propagate  together  as  a  pair.  In  some  high  mode  parameter 
regimes  there  is  an  apparent  lateral  separation  of  the  two  circulations  and  the  anti¬ 
cyclone  has  sufficient  strength  to  be  perceptible  at  the  surface.  This  is  discussed 
further  in  Appendix  I.  The  surface  cyclone  often  extends  to  depth.  In  Figure  1.8  the 
lateral  separation  is  marked  by  the  accretion  of  dye  from  the  cyclone  into  the  anti¬ 
cyclone.  There  is  frequently  some  ambiguity  in  determining  mode  number  (MIa6).  A 
half-formed  eddy  may  be  produced,  for  example,  or  there  may  be  four  distinct  waves 
present,  but  only  three  eddies  shed  due  to  coalescence.  The  measurement  uncertainty 
for  mode  number  is  then  +1. 

Our  working  definition  of  the  length  scale  of  the  eddies  is  the  average 

diameter  of  all  eddies  present.  These  measurements  were  taken  over  a  one  minute 
interval  beginning  approximately  100  s  after  the  start  of  the  experiment.  The  labo¬ 
ratory  method  of  measurement  may  be  a  consistent  underestimate.  The  eddies  must 
necessarily  be  seen  with  the  aid  of  the  source  dye.  It  is  clear  that  the  dye  cannot  be 
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transported  beyond  the  limits  of  the  circulation.  However,  it  is  not  clear  that  visible 
amounts  of  dye  reach  to  the  extent  of  the  circulation.  It  is  possible  that  the  velocity 
field  extends  beyond  the  region  of  tracer  mixing.  This  is  a  systematic  error  and  so 
will  not  affect  the  analysis  of  trends  in  the  laboratory  data,  but  it  may  be  important 
when  comparing  the  laboratory  and  numerical  data.  There  is  an  uncertainty  associ¬ 
ated  with  the  method  of  measuring  the  eddy  lengths.  This  uncertainty  comes  from 
parallax  involved  in  measuring  the  eddies  on  the  monitor  screen  instead  of  in  the 
tank,  and  with  the  resolution  of  the  picture.  The  measurements  are  good  to  t3  cm. 

The  maximum  surface  velocities  were  recorded  at  the  same  time  as 

the  eddy  sizes.  These  velocities  do  not  represent  the  maximum  velocities  which  are 
present  in  the  tank  during  the  experiment,  because  velocities  continue  to  increase 
throughout  the  duration  of  the  experiment.  In  the  laboratory  the  maximum  velocities 
were  measured  at  around  100  s  by  tracking  surface  tracers  for  3  —  45.  Sub-surface 
velocities  could  not  be  measured. 

These  speed  estimates  may  be  underestimates  for  two  reasons.  It  is  likely  that 
no  surface  tracer  was  present  where  the  velocity  had  its  maximum.  Since  a  number  of 
particles  were  tracked,  it  is  safe  to  assume  that  the  real  maximum  is  not  far  removed 
from  the  maximum  tracer  velocities,  but  it  may  be  larger.  It  is  also  possible  that  the 
tracers  themselves  introduced  some  drag,  so  that  the  velocities  measured  might  again 
be  expected  to  be  lower  than  the  actual  maximum.  The  measurement  of  velocities 
in  the  laboratory  by  tracking  the  tracers  on  the  monitor  screen  can  only  be  good  to 
within  50%  uncertainty. 

It  is  possible  to  visualize  the  shear  set  up  by  the  opposing  circulations  at  the 
surface  and  at  depth.  A  permanganate  crystal  dropped  at  the  surface  will  leave  a 
vertical  line  of  dye  as  it  descends.  In  the  region  of  the  source  the  upper  extent  of 
such  a  line  was  observed  to  move  clockwise  with  the  rotation  of  the  tank,  while  the 
lower  portion  of  the  line  is  carried  counterclockwise  away  from  it.  From  this  observed 
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shear  an  opposing  circulation  is  inferred  with  cyclonic  circulation  at  the  surface  and 
anti-cyclonic  circulation  at  depth. 

The  final  depth  of  the  convectively  produced  mixed  layer  (dl£bax)  was  recorded 
by  eye  with  a  meter  stick  held  to  the  side  of  the  tank  during  the  experiment.  Figure 
3.3  (c)  shows  a  side  view  of  convective  overturning  over  the  entire  horizontal  extent 
of  the  source.  The  deepest  point  of  the  dye  interface  was  chosen  as  dl£bax.  This  does 
not  represent  the  average  resting  depth  of  the  convectively  produced  water.  Instead, 
it  represents  the  maximum  depth  to  which  the  convection  had  an  influence.  As  will 
be  seen  below  in  4.3,  this  is  important  for  comparing  the  laboratory  and  numerical 
data.  The  measurement  uncertainty  associated  with  dl£bax  is  +  1.5  cm;  estimated  from 
the  complications  of  parallax  in  measuring  something  in  the  center  of  the  tank  from 
the  side,  and  the  ambiguity  of  the  dye  interface. 

3.3  Results 

Measurements  were  made  of  mode  number,  maximum  horizontal  speed,  eddy  length 
and  maximum  depth  in  the  laboratory  for  comparison  with  our  numerical  data  and 
to  interpret  in  the  context  of  our  scaling  laws.  In  addition,  several  more  laboratory 
experiments  were  performed  to  further  investigate  df^al.  The  laboratory  data  for 
comparison  with  the  numerical  data  is  presented  in  Table  3.2;  the  actual  comparison 
is  discussed  in  section  4.3.  Tables  3.3  and  3.4  show  the  dl£bax  data  from  the  supple¬ 
mental  experiments.  The  set  of  laboratory  and  numerical  experiments  performed  for 
comparison  have  Ba  fixed  at  .3  cm2/s3  and  R,  =  15  cm.  The  supplemental  experi¬ 
ments  were  carried  out  with  R„  =  7.5  -  22.5  cm,  to  further  examine  the  dependence 
of  depth  on  Rs. 

Figures  3.4  -  3.11  show  the  measured  values  of  Mlab,  Ul^bm,  ll“ddy  and  dl£bax.  The 
x  axis  is  the  non-dimensional  number  N/f  in  Figures  3.4,  3.6,  3.8  and  3.10,  and  R*a  in 
Figures  3.5,  3.7,  3.9  and  3.11.  The  y  axis  is  the  ratio  of  the  observed  to  the  theoretical 
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Experiment 

i 

2 

3 

4 

5 

6 

7 

8 

9 

m  i/») 

0.82 

0.45 

0.20 

0.82 

0.41 

0.22 

0.83 

0.41 

0.24 

0.50 

0.50 

0.25 

0.25 

0.25 

0.1250 

0.1250 

0.1250 

M 

3 

3 

2 

3 

2 

1 

2 

Vedduicm/s ) 

0.90 

1.40 

1.30 

1.50 

2.10 

Le<«„(cm) 

7.40 

7.60 

32.00 

29.50 

6.00 

11./16. 

23.00 

7.50 

12.25/15. /21 

27.00 

10.20 

13.50/18. 

27.00 

Table  3.2:  Laboratory  results  for  comparison  with  the  numerical  data. 


N 


0.8 

0.4 

0.2 

0.500 

4.5000 

9.0000 

17.00 

0.250 

4.00 

8.0000 

18.00 

0.125 

- 

- 

- 

Table  3.3:  Laboratory  values  for  maximum  depth,  Ra  —  7.5  cm. 


N 


0.8 

0.4 

0.2 

0.500 

- 

17.00 

- 

0.250 

- 

17.00—18.00 

- 

0.125 

- 

17.00 

- 

Table  3.4:  Laboratory  values  for  maximum  depth,  Rs  =  7.5  cm. 
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prediction.  The  uncertainty  in  measuring  the  observed  quantity  must  also  be  scaled 
by  the  same  theoretical  prediction.  Even  for  those  cases  in  which  the  measurement 
uncertainty  is  fixed,  the  predicted  value  varies  from  point  to  point,  therefore  the 
scaled  deviation  varies  as  well.  There  is  no  significant  error  associated  with  the 
theoretical  prediction.  The  figures  display  a  weighted  linear  regression  between  the 
non-dimensional  numbers  and  the  observations  normalized  by  their  predicted  values. 
Only  the  average  of  the  scaled  uncertainty  is  displayed  as  an  error  bar  in  each  figure. 
The  straight  line  displayed  is  the  fit  to  the  data.  The  dotted  lines  represent  the 
95  %  confidence  limits  associated  with  the  fit.  The  variance  of  the  fit  is  calculated 
as  a  propagation  of  the  variances  of  the  estimated  slope  and  intercept,  so  that  ajit  — 
yjc r^x2  +  <Jp.  The  dotted  lines  axe  calculated  as  the  fit  +1. 96^Jajit. 

The  non-dimensional  number  ^oes  not  introduce  any  new  parameters 

to  the  system,  therefore  it  must  be  a  linear  combination  of  the  three  non-dimensional 
numbers  with  which  we  have  chosen  to  describe  the  system.  If  the  prediction  is  an 
accurate  representation  of  the  dependence  of  observation  on  system  parameters,  then 
there  will  be  no  dependence  of  ^Action  on  any  of  the  non-dimensional  numbers.  For 
any  plot  of  vs  non-dimensional  number  there  would  be  no  trend,  ie.  a  =  0. 

The  intercept  is  the  constant  of  proportionality  between  the  observations  and  the 
predicted  scalings.  If  the  intercept  is  unity  then  all  of  the  parameters  we  are  looking 
at  have  been  accurately  portrayed  in  the  predicted  scaling. 

The  calculated  values  of  slope,  (a)  and  intercept,  (/?)  are  shown  in  Table  3.5, 
along  with  their  respective  variances  ( a2a  and  <r|).  These  values  are  plotted  in  Figures 
3.4  -  3.11.  Also  tabulated  are  the  95%  confidence  limits  for  the  slope.  Figure  3.4  shows 
plotted  against  N/f,  and  Figure  3.5  shows  versus  R*a.  The  scaling  used  here 
comes  from  the  prediction  Mscal  =  f  calculated  using  h  —  da^al  —  — • 

Recall  that  dfinal  cannot  be  measured  in  the  laboratory,  only  the  quantity  dl£ax  is 
obtained. 
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comparison 

a 

P 

a  +  e 

a  —  e 

0k  vs  N/f 

° final 

0.3605 

3.1023 

0.0053 

0.0152 

0.5033 

0.2177 

0.0137 

m.rrnmtfWiimm 

0.0041 

KiUM 

-0.0455 

mMsmm 

0.7086 

0.3028 

0.2584 

0.0033 

1.7049 

-0.2877 

KEHH 

-0.0390 

0.9900 

0.1377 

-0.2156 

-0.2900 

0.9879 

0.0708 

1.6581 

-2.2380 

mi 

-0.0948 

^||||| 

0.0032 

0.0505 

0.0159 

-0.2055 

flab 

eddy 

-1.0299 

2.6149 

0.8646 

0.1160 

Table  3.5:  Parameter  estimates  for  least  squares  fit  to  the  data.  Observational  data 
was  divided  by  scaling  prediction,  the  resulting  value  was  compared  to  the  non- 
dimensional  numbers  N/f  and  i2*.  Shown  here  are  the  estimates  of  the  slope,  a,  the 
intercept,  /3,  and  their  associated  variance,  a2a  and  ctJ,  and  the  confidence  intervals 
on  the  slope  estimate  ate  where  e  =  1.97  x  <ra,  95%  confidence.  These  quantities  are 
discussed  further  in  the  text. 

The  constant  of  proportionality  between  the  observations  and  predicted  scaling 
for  the  case  of  N/f  is  0.3204  and  for  the  case  of  R*a  is  0.3028.  As  we  will  see  in  the 
following  chapter,  the  discrepancy  between  these  constants  and  unity  can  be  explained 
by  a  missing  constant  in  the  prediction  for  df^al,  which  will  be  shown  empirically  in 
section  4.4  to  be  2.3.  When  the  missing  constant  of  proportionality  is  included, 
the  constants  are  «  0.74  and  0.70.  In  both  cases  a  zero  slope  lies  within  the  95% 
confidence  limits. 

Figure  3.6  shows  the  observations  of  maximum  horizontal  speed  against 

“rim 

N/f,  Figure  3.7  shows  the  same  quantity  against  R*0.  The  scaling  is  calculated  directly 
from  our  equation 

Kt  =  \pBM- 

The  value  ulab  =  0.9900uscai  for  the  comparison  with  N/ f  and  ulab  =  0.9546uscaf 
in  the  comparison  with  R*a.  Again  the  zero  slope  lies  within  the  confidence  limits, 
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Figure  3.4:  Laboratory  values  for  normalized  mode  number  plotted  against 

N/f. 


Figure  3.5:  Laboratory  values  for  normalized  mode  number  ( )  plotted  against 
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indicating  that  our  predicted  scaling  accurately  describes  the  dependence  of  uTim  on 
the  experimental  parameters. 

[lab 

The  next  case  is  Zeda»-  Figure  3.8  shows  the  scaled  observations  of  versus 

eddy 

N/f ,  Figure  3.9  shows  the  same  quantity  versus  R*a.  The  scaling  used  comes  from 
choosing  the  predicted  value  of  dfinai  for  h  in  our  equation 

jacal  _  Nh 
leddy  j  * 

The  value  ll^y  =  2.5842 for  the  comparison  with  N/f  and  llea^y  =  2.61491'^  in 
the  comparison  with  R*.  The  resulting  constants  of  proportionality  are  2.5842  and 
2.6149.  The  discrepancy  between  these  and  unity  can  be  completely  explained  again 
by  the  factor  of  2.3  associated  with  dfinal.  The  corrected  values  are  «  1.123  and 
1.136.  In  both  cases  zero  slope  lies  within  the  confidence  limits. 

The  last  case  is  dmax.  Figure  3.10  shows  Jlfff-  versus  N/f,  while  Figure  3.11 

/trial 

shows  the  same  quantity  versus  R*a.  In  neither  case  does  the  zero  slope  lie  within 
the  confidence  limits.  This  suggests  that  the  laboratory  measurement  of  dl£fax  is  not 
accurately  described  by  the  prediction  for  dsf:ffal. 
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Figure  3.7:  Laboratory  values  for  normalized  maximum  velocity  plotted  against 


Figure  3.8:  Laboratory  values  for  normalized  eddy  length  plotted  against  N/f. 

eddy 
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eddy  length  to  predcted  scaling 


Figure  3.9:  Laboratory  values  for  normalized  eddy  length  (n^)  plotted  against  R*. 

**4dy 
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Figure  3.11:  Laboratory  values  for  normalized  maximum 

K 


depth  (-jffff  )  plotted  against 

final 
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Chapter  4 


The  Computational  Analog 

4.1  Numerical  parameters 

The  numerical  model  used  in  this  study  is  the  non-hydrostatic  primitive  equation 
model  described  in  detail  in  [Brugge  et  al.,  1991].  The  model  was  developed  from  the 
atmospheric  model  written  by  Miller  (1971),  and  adapted  to  become  a  model  appro¬ 
priate  for  studying  ocean  convection.  The  equations  for  a  Boussinesq  incompressible 
fluid  with  an  equation  of  state  appropriate  for  sea-water  are 

Momentum: 

( d2u  d2u\  d2u 

VH  \d^+d^J  +Vvd * 

(d2v  d2v\  d2v 

VH  \d^  +  d^J  +Vvdz* 

( d2w  d2w\  d2w 

VH  (a?+ 


Du 

+ 

1  dp' 

fv 

d] 

Po  dx 

Dv 

7>i 

+ 

1  dp' 
Po  dy 

+ 

fu 

Dw 

Dt 

+ 

1  dp' 
po  dz 

+ 

p' 

9 

Po 
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Continuity: 


Density  Conservation: 


du  dv  dw 
dx  dy  ^  dz 


Equation  of  State: 


P  =  Pail  -  -  s„)) 

Where: 

D  d  d  d  d 

— -  =  — +  «- — |-  v— — \-w-p—  m  height  coordinates. 

Dt  at  ox  dy  dz 

Here  u,  v  and  w  are  the  customary  velocities  in  the  x,  y  and  z  directions,  Q„ 
is  equivalent  salt  flux,  p  is  density,  p  is  pressure,  and  vh  and  i/y  are  the  horizontal 
and  vertical  coefficients  of  diffusion  respectively. 

The  full  equation  of  state  p  =  p  (T,S,z)  is  not  required  in  the  parameter  regime 
of  the  laboratory  as  pressure  will  never  be  significant  in  the  range  of  30-40  cm  depth, 
and  temperature  never  varies  from  21  °C. 

In  these  equations  p'  is  the  perturbation  of  the  density  field  from  its  initial 
(reference)  state  ps ;  p0  is  a  standard  constant  value.  Similarly,  p’  is  the  deviation  of 
pressure  from  the  initial  hydrostatically  balanced  p3. 

The  model  is  doubly  periodic  and  has  a  rigid  lid  and  free-slip  bottom.  For 
practical  purposes  this  means  that  there  is  no  side  wall  interference  until  a  migrating 
eddy  leaves  the  domain  to  one  side  and  re-enters  from  the  other  side.  There  will  be  no 
bottom  drag  to  consider.  The  rigid  lid  also  ensures  that  there  are  no  surface  gravity 
wave  modes  to  cause  numerical  problems. 

A  time  step  of  0.1  s  is  used  to  satisfy  the  CFL  conditions  with  grid  spacing  of 
only  0.86  and  1.56  cm  in  the  horizontal  and  vertical  respectively. 


44 


The  buoyancy  forcing  is  applied  over  a  circular  patch  of  desired  radius.  Effec¬ 
tively  the  model  introduces  salt  to  the  appropriate  disk  of  water,  one  grid  space  deep. 
There  is  no  addition  of  volume  to  the  system,  only  salt. 

The  physical  values  of  kinematic  diffusivity  were  used  in  the  model  v  = 
10-6  m2/s.  The  molecular  diffusion  coefficients  for  salt  in  water  were  adjusted  in 
keeping  with  stability  criteria.  Instead  of  the  molecular  value  k„  =  10-9  m2/s,  the 
model  was  initiated  with  kv  —  Kh  —  10-6  m2/s,  where  kv  and  Kh  are  the  vertical  and 
horizontal  diffusion  coefficients  respectively. 


4.2  Measurement  and  qualitative  results 

The  numerical  experiments  were  set  up  to  first  duplicate  the  laboratory  experiments 
and  then  expand  the  ranges  of  N ,  /,  and  Rs  beyond  the  capabilities  of  the  laboratory. 
Table  4.1  lists  the  experimental  parameters  for  each  run.  Also  displayed  are  the  non- 
dimensional  numbers  N/ f,  A*,  and  Rs/H,  where  H  is  again  30  cm  and  is  the  total 
depth  of  the  numerical  tank.  The  first  set  of  runs  was  performed  with  a  constant 
value  of  R3  =  15  cm,  values  of  N  ranging  from  0.18  <  N  <  0.9  /s  and  values  of  f 
ranging  from  0.125  <  /  <  1.0  /s.  The  second  set  of  experiments  were  performed 
varying  the  value  of  Ra  from  7.5  <  Rs  ^  22.5  cm,  while  keeping  constant  values 
of  N  «  0.4  /s  and  /  =  0.25  /s. 

The  values  measured  in  the  numerical  work  are  mode  number  of  the  baroclinic 
instability  ( Mnum ),  length  scale  of  the  baroclinic  eddies  maximum  velocity 

final  resting  depth  of  the  convectively  produced  water  ( ) ,  maximum 
depth  of  convective  influence  for  comparison  with  the  laboratory  and  final 

volume  of  convectively  produced  water  (Vp£™i). 

Mode  number  (Mnum)  is  measured  in  the  computations  the  same  way  it  is 
measured  in  the  laboratory.  However,  instead  of  the  dye  field,  the  velocity  field  is 
used.  For  future  experiments  it  would  be  possible  to  advect  a  passive  tracer  in  the 
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Table  4.1:  This  table  displays  the  experimental  parameters  for  the  numerical  exper¬ 
iments.  Also  shown  are  the  values  of  the  significant  non-dimensional  parameters. 
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numerical  experiments.  Figure  4.1  (a)  shows  an  unobscured  overhead  view  of  the 
velocity  field  at  the  beginning  of  an  experiment,  after  the  forcing  has  been  initiated. 
The  plot  shown  is  a  horizontal  cross  section  at  level  2  in  the  model,  the  level  at  the 
bottom  of  the  source  block.  Figure  4.1  (b)  shows  a  side  view  of  the  experiment;  the 
picture  displays  a  cross  section  through  the  center  of  a  chimney  where  the  contours 
depict  salt.  Figure  4.1  (c)  shows  a  typical  mode  three  instability,  displaying  several 
well  developed  eddies  propagating  away  from  the  source.  Figure  4.1  (d)  shows  a 
high  mode  instability,  illustrating  the  surface  manifestation  of  the  sub-surface  anti¬ 
cyclone.  The  anti-cyclonic  signature  is  discernible  because  of  its  strength  and  the 
lateral  separation  of  the  two  circulations. 

The  length  scale  of  the  eddies  (l™™)  is  operationally  defined  as  the  average 
diameter  of  all  eddies  present.  These  measurements  were  made  for  comparison  with 
the  laboratory  data.  In  the  numerical  experiment  the  eddies  were  measured  at  one 
time  only,  100  s  after  the  start  of  the  experiment.  We  have  discussed  above  the 
fact  that  the  eddies  continue  to  grow  in  time.  Therefore  the  only  way  to  accurately 
compare  eddies  from  two  different  experiments  is  to  measure  them  at  the  same  time. 
The  edge  of  an  eddy  was  marked  by  the  abrupt  decrease  in  the  velocity  field.  As 
discussed  above,  it  is  possible  that  the  velocity  field  extends  beyond  the  region  of 
tracer  mixing,  potentially  making  the  numerical  observations  systematically  larger 
than  the  laboratory  observations. 

The  maximum  surface  velocities  (u?™) were  recorded  at  the  same  time  as 
the  eddy  sizes,  for  comparison  with  the  laboratory  data.  As  mentioned  above,  the 
laboratory  measurements  may  be  an  underestimate  of  the  maximum  velocity  present. 
In  the  numerical  model  it  is  possible  to  locate  the  grid  point  with  the  maximum 
velocity  present  in  the  field  at  a  given  time. 

The  final  depth  of  the  convectively  produced  water  was  recorded  using 

information  from  volume  histograms.  The  time  evolution  of  the  volume  of  water 
within  a  density  range  can  give  important  insights  into  the  dynamics  of  the  convective 
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Figure  4.1:  Horizontal  and  vertical  cross  sections  taken  from  the  numerical  model: 
a)  overhead  view  of  numerical  experiment  at  start,  note  the  formation  of  individual 
overturning  plumes;  b)  side  view  of  the  numerical  experiment,  note  the  narrow  chim¬ 
ney  beneath  the  source,  and  the  spreading  of  the  dense  fluid  at  depth;  c)  overhead 
view  of  numerical  experiment,  distinct  eddies  have  formed  and  are  propagating  away 
from  the  source;  d)  overhead  view  at  surface;  in  this  case,  the  anti-cyclonic  circulation 
is  apparent  at  the  surface. 
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event.  As  discussed  above,  the  laboratory  measurement  of  dl^ax  gives  information  only 
on  the  deepest  point  that  had  contact  with  the  convection  event  at  any  time.  For  the 
scaling  arguments  it  is  preferable  to  know  the  average  resting  depth  of  the  convectively 
produced  water.  For  the  sake  of  comparison  it  is  necessary  also  to  measure  the 
numerical  analog  to  the  laboratory  observation,  a  property  here  assigned  as 
The  volume  histograms  are  used  to  obtain  values  of  dj^u  Ci  and  ■ 

Figure  4.2  illustrates  the  volume  histogram  time  series.  The  fraction  of  the 
total  volume  lying  within  a  given  density  range  can  be  calculated  from  the  salt  field  of 
the  model.  The  value  of  salinity  is  known  at  every  grid  point  for  every  time.  If  density 
is  binned  to  reflect  the  initial  stratification,  then  a  histogram  at  t=0  s  will  show  a 
flat  line,  indicating  an  equal  amount  of  water  in  each  density  bin.  This  corresponds 
to  the  initialization  of  each  layer  in  the  absence  of  forcing.  The  flat  line  in  Figure  4.2 
shows  5.26%  of  the  total  volume  of  water  in  each  of  the  19  bins,  corresponding  to  the 
19  levels  used  in  the  model. 

At  t — 0. 1  s  the  buoyancy  forcing  is  switched  on.  Then  surface  water  is  made 
more  dense.  As  mentioned  above,  the  buoyancy  forcing  in  the  model  takes  the  form 
of  added  salinity,  with  no  increase  in  the  volume  of  water  in  the  domain.  The  surface 
water  is  converted  to  denser  water.  This  is  manifest  in  the  first  peak  in  Figure  4.2. 
This  line  corresponds  to  the  histogram  at  t=25  s.  The  amount  of  water  in  bins  1 
and  2  has  decreased,  illustrating  the  conversion  of  surface  water  to  denser  water. 
This  water  shows  up  as  a  surplus  in  bins  3,4  and  5.  The  next  peak  corresponds  to 
t=50  s.  There  is  a  larger  surplus,  now  in  bins  5,6  and  7.  More  surface  water  has 
been  converted,  and  all  of  the  mixed  water  has  been  subjected  to  more  buoyancy 
forcing.  The  result  is  a  larger  mass  of  water  at  a  higher  density.  As  buoyancy  forcing 
continues,  this  volume  continues  to  increase.  Note  however,  that  the  density  of  the 
convectively  produced  water  does  not  continue  to  increase  indefinitely.  In  every  case 
studied  the  final  density  was  achieved  within  the  first  300  s. 
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Figure  4.2:  Shown  here  is  a  series  of  volume  histograms  for  a  typical  numerical 
experiment.  Each  curve  represents  a  histogram  at  a  particular  time,  during  t=0- 
300  s.  The  flat  line  indicates  the  t=0  s  initial  stratification  with  no  modification  to 
any  water  mass.  The  first  peak  represents  t=25  s.  The  intermediate  peak  is  from 
t=50  s.  The  last  peak  is  from  t=300  s.  The  estimates  for  Vfinai,  djinai  and  dmax  were 
all  measured  at  this  time. 
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Experiment 

1 

2 

3 

4 

5 

6 

am 

8 

9 

N(l/s) 

0.9370 

0.4310 

0.24 

0.9370 

0.4310 

0.24 

0.9370 

0.4310 

m/s) 

0.50 

0.50 

EH 

0.25 

0.1250 

0.1250 

M 

5 

6 

4 

4 

4 

3 

2 

3 

2 

Veddy{cm/s) 

1.71 

2.48 

2.14 

2.41 

2.79 

2.83 

3.00 

2.83 

Lcddy{cm) 

15.70 

14.80 

- 

11.87 

- 

- 

- 

7.89 

14.21 

EH 

18.95 

7.89 

14.21 

17.37 

Table  4.2:  Numerical  results  for  comparison  with  the  laboratory  data. 

Now,  several  useful  quantities  can  be  defined.  The  depth  at  which  water  of 
that  density  would  be  neutrally  buoyant  in  the  background  stratification  is  defined 
as  the  final  depth  of  that  water  The  density  of  the  deepest  zero  crossing  of 

the  convectively  produced  water  is  defined  as  the  maximum  density  effected  by  the 
mode  water.  The  depth  at  which  water  of  that  density  would  be  neutrally  buoyant  is 
defined  as  the  maximum  depth  of  influence  of  the  convection  (d™™).  The  final  volume 
of  convectively  produced  water  ( v "JJJJjj)  is  defined  to  be  any  water  occupying  the  peak 
densities  in  excess  of  the  original  stratification.  That  is,  the  original  stratification 
is  subtracted,  and  an  integral  is  performed  over  the  dominant  peak,  as  indicated  in 
Figure  4.2.  All  of  these  quantities  were  measured  at  t  =  300  s. 


4.3  Comparison  with  the  laboratory 

As  stated  in  section  3.3,  laboratory  measurements  were  made  of  mode  number,  max¬ 
imum  horizontal  speed,  eddy  length  and  maximum  depth,  for  comparison  with  the 
numerical  data.  The  numerical  estimates  of  these  properties  are  presented  in  Table 
4.2.  Figure  4.3  shows  the  observed  mode  number,  M  from  the  laboratory  plotted 
against  the  numerical  estimates.  The  error  bars  in  x  and  y  indicate  the  measurement 
uncertainty  associated  with  the  numerical  and  laboratory  estimates  respectively.  All 
of  the  data  lie  within  an  error  bar  of  a  one  to  one  relationship.  This  indicates  that 
the  numerics  are  properly  resolving  the  baroclinic  instability  wavelength. 
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Laboratory  estimate  of  mode  number 


Numerical  estimate  of  mode  number 


Figure  4.3:  Comparison  between  laboratory  and  numerical  estimates  of  mode  number, 
Mlab  and  Mnum . 
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Numerical  estimate  of  velocity 


Figure  4.4:  Comparison  between  laboratory  and  numerical  estimates  of  maximum 
velocity,  u1^  and  u™™. 

Figure  4.4  shows  the  observations  of  maximum  velocity,  from  the  laboratory 
urirn  versus  numerics  The  estimates  from  the  numerical  experiments  are  consis¬ 

tently  higher  than  the  estimates  from  the  laboratory.  This  eventuality  was  discussed 
above  as  a  possible  consequence  of  the  laboratory  measurement  errors.  The  fact  that 
the  data  are  so  closely  grouped  and  the  consistency  of  the  small  discrepancy  seem  to 
indicate  a  satisfactory  agreement  between  the  laboratory  and  numerical  estimates  of 
horizontal  speed. 

Figure  4.5  shows  the  laboratory  and  numerical  observations  of  eddy  length 
scale  ll“£dy  and  Again  as  expected,  the  numerical  estimate  is  consistently  larger 
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than  the  laboratory  estimate.  This  discrepancy  is  much  larger  than  that  in  the  ve¬ 
locity  data.  The  discrepancy  is  10  —  15cm.  In  light  of  the  good  linear  relationship 
between  the  laboratory  and  numerical  estimates  it  seems  more  likely  to  be  another 
symptom  of  systematic  measurement  error  rather  than  a  symptom  of  numerical  fail¬ 
ure.  As  discussed  in  section  3.2,  the  laboratory  eddies  can  only  be  measured  to  the 
extent  that  a  visible  dye  signature  exists.  This  dye  edge  may  lie  closer  to  the  center 
of  the  eddy  than  the  edge  of  the  azimuthal  velocity  field  measured  in  the  numerical 
experiment.  Any  lag  in  the  time  of  measurement  of  the  eddies  between  the  labora¬ 
tory  and  numerics  would  only  aggravate  the  issue.  The  close  grouping  of  the  data 
about  a  linear  relation  indicates  a  satisfactory  agreement  between  the  laboratory  and 
numerical  data. 

Finally,  Figure  4.6  shows  the  comparison  between  the  laboratory  and  numerical 
estimates  of  dl^ax  and  d™™.  This  figure  4.6  shows  that  the  agreement  between  the 
two  estimates  is  approximately  one  to  one,  within  measurement  error.  There  is  some 
discrepancy  at  large  depths.  This  may  come  from  the  systematic  correlation  between 
the  numerical  estimates  of  and  We  have  already  established  that  the 

laboratory  estimate  of  d!^ax  does  not  scale  with  ,  but  the  numerical  estimate 
of  does  better,  as  we  will  see  in  the  following  section.  So,  the  discrepancy  is 
not  surprising,  it  indicates  that  the  numerical  estimate  of  d%f™  is  not  a  perfect  proxy 
for  the  laboratory  estimate  d!^bax.  However,  the  discrepancy  is  small  enough  for  us 
to  draw  some  useful  information  from  the  comparison.  The  adequate  agreement  is 
significant  in  that  it  supports  the  supposition  that  the  individual  plumes  need  not  be 
perfectly  resolved  in  order  to  accurately  simulate  chimney  evolution.  The  laboratory 
presents  us  with  a  physical  system  that  must  evolve  within  any  existing  influence  of 
the  plumes.  The  numerical  experiments  were  not  adequately  resolving  the  smallest 
of  the  individual  plumes.  Yet,  the  vertical  evolution  of  the  chimney  is  in  each  case 
virtually  identical  within  measurement  error. 
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Laboratory  estimate  of  maximum  depth 


Numerical  estimate  of  maximum  depth 


Figure  4.6:  Comparison  between  laboratory  and  numerical  estimates,  dl£bax  and 
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4.4  Results 


In  addition  to  the  measurements  made  for  comparison,  numerical  estimates  were  made 
of  over  a  range  of  Rs  values  and  a  wider  range  of  N  and  /.  At  the  same  time, 
estimates  were  made  of  and  the  final  volume  of  dense  water  produced 

by  the  convective  event.  These  quantities  are  presented  in  Table  4.3,  along  with 
the  experimental  parameters  which  were  varied;  the  table  also  repeats  the  values 
of  used  above  for  comparison  with  the  laboratory.  For  all  these  experiments 

Ba  —  .3cm2 /s3,  v  —  10 ~6m2/s  and  k  =  10 ~6m2/s. 

Again,  as  in  section  3.3,  the  observations  are  scaled  by  prediction  and  presented 
against  the  non-dimensional  numbers.  Now  we  have  a  range  of  Ra/H  to  explore  also. 
The  error  analysis  is  identical  to  that  in  Chapter  3,  with  the  exception  of  the  Vfinai 
analysis  which  will  be  discussed  in  detail  below.  The  values  of  slope,  (a)  and  intercept, 
(/3)  are  shown  in  Table  4.4,  along  with  their  respective  variances  (a2  and  cr|).  Also 
tabulated  are  the  95%  confidence  limits  for  the  slope. 

jnum 

Figures  4.7,  4.8,  and  4.9  show  the  scaled  quantity  ^ r  against  AT//,  R*,  and 

a  final 

Raj H  respectively.  The  scaling  used  here  is 

a  jy 

In  all  three  cases  the  zero  slope  lies  within  the  confidence  limits. 

The  resulting  constants  of  proportionality  between  the  observed  and  predicted 
values  are  2.2,  2.3  and  2.1.  The  arguments  which  led  to  the  relation  for  were  not 
capable  of  determining  this  constant.  The  reader  should  recall  that  our  derivation 
of  this  quantity  came  from  an  assumption  that  all  of  the  potential  energy  in  the 
baroclinic  zone  would  be  converted  to  kinetic  energy.  Our  empirical  constant  is  an 
indicator  of  the  efficiency  of  that  conversion. 

Figures  4.10,  4.11  and  4.12  show  the  quantity  Jiff-  versus  N/ /,  i?*,  and  Rsj  H. 

final 

The  zero  slope  lies  within  the  confidence  limits  only  for  the  cases  of  R*  and  Rs/H. 
In  the  case  of  N/f  there  exists  a  significantly  non-zero  positive  relationship.  In  light 
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Comp  run 

f 

N 

R,  (cm) 

^macc(c^) 

d/ina/(c7Tl) 

Volftnal(%) 

1 

0.5000 

0.9370 

7.89 

4.74 

1.08 

2 

m 

11.05 

7.89 

1.62 

3 

0.2400 

7.5 

15.79 

2.46 

4 

BKTO1 

7.5 

6.32 

1.43 

_ 

5 

0.2500 

7.5 

11.05 

7.89 

6 

0.2500 

n 

7.5 

14.21 

2.67 

7 

0.1250 

1  .  : 

™ _ 

| 

0.87 

8 

0.1250 

0.4310 

9.47 

1.33 

9 

0.1250 

0.2400 

7.5 

14.21 

2.62 

10 

0.5000 

7.89 

2.34 

11 

0.5000 

WWW 

14.21 

5.15 

12 

0.5000 

15.0 

18.95 

6.25 

13 

0.2500 

0.9370 

15.0 

7.89 

6.32 

2.86 

14 

mm 

14.21 

4.67 

15 

0.2400 

EH 

18.95 

7.75 

16 

0.9370 

7.89 

4.74 

4.52 

17 

14.21 

11.05 

4.84 

18 

17.37 

14.21 

8.95 

19 

0.2500 

0.1800 

20.53 

17.37 

9.72 

20 

1.0000 

0.2400 

20.53 

7.68 

21 

0.2500 

0.4310 

7.62 

22 

0.5000 

0.4310 

22.5 

11.18 

23 

0.1250 

22.5 

9.47 

4.45 

Table  4.3:  This  table  displays  the  raw  data  for  the  numerical  experiments  along  with 
the  experimental  parameters  used  for  each  run. 
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arison 

a 

0 

oa 

Op 

a  +  e  a  —  e 

vs  N/f 

0.0882 

2.2025 

0.0040 

0.0116 

0.2127  -0.0362 

vs  Rl 

-0.3614 

2.3807 

0.2395 

0.0123 

0.5978  -1.3206 

vs  R,/H 

0.3825 

2.1379 

0.2296 

0.0551 

1.3216  -0.5566 

maim 

0.2882 

2.5550 

0.0040 

0.0116 

0.4127  0.1638 

VS  R*0 

-0.1222 

2.9502 

0.2395 

0.0123 

0.8369  -1.0814 

vs  R,/H 

0.2266 

2.8226 

0.2296 

0.0551 

1.1658  -0.7125 

Table  4.4:  Measured  parameters  of  least  squares  fit  to  the  data.  Observational  data 
was  divided  by  scaling  prediction,  the  resulting  value  was  compared  to  the  non- 
dimensional  numbers  N/f  and  R*a.  Shown  here  are  the  estimates  of  the  slope,  a,  the 
intercept,  j3,  and  their  associated  variance,  oa  and  op,  and  the  confidence  intervals 
on  the  slope  estimate  ale  where  e  =  1.97  x  o]J2. 
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Figure  4.8:  Numerical  value  of  .ffi1  plotted  against  R*. 
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Figure  4.9:  Numerical  value  of  plotted  against  RJH. 
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of  the  successful  agreement  of  our  estimate  of  with  scaling,  the  discrepancy 
between  and  d'f^al  seems  to  indicate  an  inadequacy  in  our  definition  of  d%/™,  or 
our  choice  of  as  a  proxy  for  dJgE,,  as  opposed  to  an  inadequacy  in  the  scaling. 
Prom  Table  4.3  we  see  that  the  constants  of  proportionality  relating  to  dsfffal 
are  slightly  larger  than  those  relating  with  djw,,  again  owing  to  definition. 

Figures  4.13,  4.14  and  4.15  show  the  measured  quantity  Vf divided  by  the 
predicted  value  of  Vfgj  =  TrR4J3dfinaiB^3  x  t,  versus  N/f,  R*0,  and  R,/H.  As  stated 
above,  the  values  of  were  measured  only  at  one  time;  therefore  the  appropriate 
scaling  depends  on  a  constant  value  of  time,  incorporated  here  in  the  constant  of 
proportionality.  This  fit  tests  only  the  dependence  of  on  the  predicted  rate  of 
production.  There  is  no  significant  measurement  uncertainty  associated  with  estimate 
of  volume.  A  non-weighted  linear  regression  on  this  data  indicates  that  there  is  no 
significant  non-zero  trend  in  any  of  the  three  cases,  with  95  %  confidence.  In  the 
comparison  with  N/f  the  estimated  slope  is  a  =  1.7846  and  the  estimated  intercept 
is  P=  19.0706. 
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Figure  4.11:  Numerical  value  of  Jiff-  plotted  against  R*. 
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Figure  4.12:  Numerical  value  of  plotted  against  RJH. 
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Figure  4.13:  Numerical  value  of  plotted  against  AT//. 
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In  the  comparison  with  i?*  the  estimated  slope  is  a  —  6.3557  and  the  estimated 
intercept  is  /?  =  21.8819. 

In  the  comparison  with  R„/H  the  estimated  slope  is  a  =  —18.2357  and  the 
estimated  intercept  is  (3  =  31.0357.  Again,  in  each  of  the  three  cases,  the  slope  is  not 
statistically  different  from  zero  within  95  %  confidence. 
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Figure  4.14:  Numerical  value  of  plotted  against  jR*. 


Figure  4.15:  Numerical  value  of  plotted  against  Ra/H 


Chapter  5 


Discussion 


The  quantities  estimated  from  the  numerical  model  compared  satisfactorily  with  the 
laboratory  observations.  There  were  some  discrepancies,  but  none  which  contradicted 
the  premise  that  the  numerical  model  is  adequately  resolving  the  relevant  physics  we 
have  chosen  to  consider  in  our  problem. 

The  laboratory  and  numerical  experiments  agree  best  in  estimates  of  mode 
number.  The  laboratory  observations  of  this  quantity  agree  very  well  with  the  pre¬ 
dicted  scaling.  This  was  expected,  given  that  mode  number  is  well  understood.  How¬ 
ever  it  is  encouraging  in  its  validation  of  the  numerical  model. 

The  estimates  of  maximum  velocity  from  the  laboratory  and  numerics  agree 
well  with  each  other,  although  the  laboratory  measurement  is  a  consistent  under¬ 
estimate  of  the  numerical  value.  The  laboratory  observations  are  associated  with 
very  large  measurement  uncertainty.  Within  the  limits  of  that  uncertainty  they  agree 
excellently  with  the  predicted  scaling. 

The  numerical  and  laboratory  estimates  of  leddy  also  agree  well,  while  the 
laboratory  definition  is  again  an  underestimate  of  the  numerical  value.  The  laboratory 
estimates  of  agree  well  with  the  predicted  scaling  only  if  one  chooses  h  ~  dj^al 
as  the  appropriate  scale  for  depth.  That  is,  one  must  use  the  calculated  prediction  of 
dSfinai  since  there  is  no  laboratory  observation  of  this  quantity.  The  measurement  of 
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is  n°t  an  adequate  proxy.  If  one  uses  h  ~  dl£bax  then  ll^dy  does  not  agree  with 
the  predicted  scaling. 

The  laboratory  and  numerical  estimates  of  dmax  did  not  compare  as  well  as  the 
estimates  of  other  quantities.  This  seems  to  be  an  artifact  of  the  definitions  of 
and  dTfiwi-  The  laboratory  estimates  of  dl£bax  did  not  agree  with  the  predicted  scaling 
for  dy^ai ■  The  laboratory  measure  of  dl£bax  reflects  all  overshooting  and  oscillation 
past  the  spreading  level.  This  is  not  the  quantity  we  actually  wanted  to  measure. 
We  are  interested  in  the  level  of  the  bulk  of  the  dense  water.  To  observe  this  more 
accurately,  the  laboratory  apparatus  could  have  been  left  running  for  several  hours 
after  the  forcing  was  terminated.  The  convectively  produced  water  would  gradually 
rebound  to  its  neutrally  buoyancy  level.  A  measurement  at  that  time  would  more 
closely  estimate  our  quantity  d‘^^[al . 

The  numerical  estimates  of  agree  well  with  the  predicted  scaling  for  that 
quantity.  This  agreement  supports  our  choice  of  L  =  as  the  length  scale  for  Ay 
the  distance  across  the  baroclinic  zone.  It  also  validates  our  assumption  of  thermal 
wind  balance  in  the  convective  region.  It  is  interesting  to  note  that  while  we  are 
describing  the  impact  of  baroclinic  instability  on  mixed  layer  depth,  the  relation  does 
not  depend  explicitly  on  f.  According  to  our  results,  the  final  depth  of  the  convec¬ 
tively  produced  water  depends  only  on  the  size  and  intensity  of  the  forcing  and  the 
ambient  stratification.  Implicit  in  the  derivation  of  this  quantity  da^al  is  the  balance 
between  buoyancy  forcing  at  the  surface  and  the  lateral  transport  of  dense  fluid  by 
baroclinic  eddies.  Again  recall  that  one  of  the  most  important  ramifications  of  the 
calculated  constant  of  proportionality  for  is  an  indication  of  the  efficiency  of 

the  baroclinic  transport.  The  quantity  was  derived  from  the  energy  balance  assuming 
that  all  available  potential  energy  was  converted  to  kinetic  energy.  The  constant  of 
proportionality  between  df“™[  and  d3J^al  is  a  gauge  of  the  accuracy  of  that  assump¬ 
tion.  This  balance  then  is  important  not  only  in  the  matter  of  convection,  but  also 
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anywhere  that  the  efficiency  of  baroclinic  eddies  is  of  interest,  for  example  in  the 
parameterization  of  baroclinic  eddies  for  use  in  coarse-scale  models. 

In  the  case  of  convection,  where  forcing  persists  long  enough  for  the  forma¬ 
tion  of  baroclinic  eddies,  the  scales  summarized  here  will  govern  the  evolution  of  a 
convective  chimney: 

Macai  — 

Nd. final 5 

jscal  _  Ndfjnat 
* eddy  y  ) 

and 

a  final  jy 

The  data  from  both  laboratory  and  numerical  experiments  support  these  scal¬ 
ings.  To  see  when  these  scalings  would  be  relevant  we  need  to  take  a  closer  look  at  the 
magnitude  of  “long  enough” .  Recall  that  the  shortest  time  necessary  for  the  chimney 
to  evolve  to  dfinai  is  given  by 


t 


M  ( 1/3 

\/2  VBoJ 


using  our  empirical  constant  of  proportionality  for  dfinai  and  the  y/2  from  the  one 
dimensional  derivation.  Using  our  above  values  for  a  typical  Meditteranean  event 
B0  «  800 W/m2  over  an  area  Ra  —  20  —  45 km  we  find  t  =  2  x  105s,  or  «  5days.  The 
Mistral  typically  persist  for  up  to  a  week  however,  forcing  this  strong  rarely  persists 
for  more  than  1-2  days.  Therefore,  it  is  unlikely  that  baroclinic  eddies  are  impacting 
the  evolution  of  Meditteranean  deep  water. 

Now  let  us  look  at  some  typical  values  for  the  Labrador  Basin.  We  will  consider 
the  conditions  at  weather  station  Bravo  as  typical  of  the  entire  basin.  There  is 
some  evidence  that  the  buoyancy  flux  measured  here  is  typical  of  the  basin  based 
on  analyzed  fields  from  European  Center  for  Medium-Range  Weather  Forecasting 
(ECMWF),  [Pickart,  1994].  Then  choosing  average  winter  values  from  the  Bravo 
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records  for  several  years  we  find  for  warm  winters  a  basin  wide  heat  loss  of  100 W/m2 , 
and  for  extreme  cold  winters  values  of  600 Wjm2.  These  heat  losses  persist  for  a 
month  or  more. 

Using  a  value  of  buoyancy  frequency  calculated  from  CTD  data  [Pickart,  1994] 
it  is  reasonable  to  assert  that  the  ambient  stratification  encountered  by  the  convective 
overturning  is  N  «  9  x  10_4s_1.  Using  these  values  and  the  dimensions  of  the  basin 
(*.  500&m)  we  find  the  case  of  B0  =  100 W/m2,  to  corresponds  with  a  depth 

d final  =  690m  with  a  minimum  formation  time  of  t finai  =  33.5 days.  The  case  of 
B0  =  600 W/m2  corresponds  with  a  final  depth  dfinai  =  1260m  and  a  minimum 
formation  time  of  tfinai  —  18.5 days.  Note  that  the  stronger  the  forcing  the  faster  the 
convective  penetration  and  the  deeper  the  mixing. 

The  depths  predicted  above  are  roughly  consistent  with  recorded  depths  of 
mixing  at  Bravo  [Lazier,  1980].  For  example,  the  estimate  of  100 W/m2  is  appropriate 
to  winter  1968,  where  the  recorded  depth  of  mixing  was  800m  and  we  predicted  690m. 
Winter  1972  was  extremely  cold,  Ba  «  800 W/m2.  That  winters  recorded  depth  of 
mixing  was  1500m  where  we  predicted  1260m. 

Given  these  estimates  it  seems  likely  that  Labrador  Sea  overturning  is  affected 
by  the  processes  studied  here.  In  that  case,  baroclinic  eddies  formed  at  the  edges  of 
the  convection  site,  ie.  the  edge  of  the  basin,  provide  an  appealing  mechanism  for  the 
transfer  of  convectively  produced  water  into  the  DWBC. 

One  might  also  expect  that  the  final  volume  of  convectively  produced  water 
manufactured  by  one  of  these  events  would  also  be  governed  by  the  above  scaling 
theories.  The  numerical  estimates  of  Vfinai  agreed  well  with  the  proposed  dependence 
on  rate 

D5/3  d2/3 

ratemax  =  R4/3  B^d^i  =  *  N°  , 

recall  that  the  final  volume  is  assumed  to  depend  on  the  maximum  rate  of  production 
and  the  duration  of  the  forcing 

Vfinai  =  TQ,t(imax  * 
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where  our  measurements  were  taken  at  t  —  300  s.  This  scaling  indicates  that  the 
final  volume  of  water  produced  will  depend  on  the  size  of  the  forcing  area  and  the 
final  depth  of  penetration,  or  more  explicitly  it  will  depend  on  the  size  and  intensity 
of  the  forcing,  and  the  initial  stratification. 

The  final  volume  of  water  produced  is  of  interest  in  observations  of  the  Labrador 
Sea.  There  exist  volume  estimates  of  water  produced  by  convection  in  the  Labrador 
Sea.  This  Labrador  Sea  Water  (LSW)  is  associated  with  the  Deep  Western  Boundary 
Current  (DWBC)  as  it  leaves  the  Labrador  Basin,  and  transports  have  been  calcu¬ 
lated  from  hydrographic  sections  through  the  current.  However,  very  little  is  known 
about  the  frequency,  duration  or  position  of  convective  events  in  the  Labrador  Sea. 
It  is  not  clear  that  water  with  the  characteristics  of  LSW  is  produced  every  year; 
but  the  signal  of  LSW  associated  with  the  DWBC  is  relatively  steady.  In  unraveling 
these  questions  it  will  be  useful  to  understand  how  much  water  to  expect  from  a 
given  convection  event;  and  to  know  that  the  estimate  may  be  made  given  a  hydro- 
graphic  survey  of  the  ambient  stratification  of  the  region,  and  meteorological  surveys 
of  typical  storms. 
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Appendix  I 

There  is  one  other  point  of  interest  regarding  the  occurrence  of  heton  pairs  in  our 
experiments.  The  term  heton  pair  refers  to  the  phenomenon  of  lateral  separation  of 
the  surface  cyclone  and  the  subsurface  anti-cyclone.  In  some  regimes  we  see  a  baro- 
clinic  eddy  whose  sole  surface  signature  is  cyclonic,  with  an  anti-cyclonic  circulation 
present  only  at  depth.  Such  cases  can  be  seen  in  several  of  our  numerical  and  labo¬ 
ratory  experiments  and  in  oceanic  regimes  described  in  [Jones  and  Marshall,  1993]. 

However,  in  regimes  of  high  mode  instability  we  see  the  lateral  separation  of 
the  two  circulations,  and  a  surface  signature  associated  with  the  anti-cyclone.  The 
cyclonic  circulation  is  always  the  strongest  at  the  surface,  and  is  often  the  stronger 
of  the  pair,  extending  to  depth.  The  anti-cyclone  is  always  strongest  at  depth  and  is 
occasionally  too  weak  to  be  visible  at  the  surface.  Figure  5.1  displays  a  map  of  the 
numerical  experiments  in  which  the  lateral  separation  of  the  eddies  was  manifest. 


74 


Presence  of  heton  pairs  in  the  N— f  plane 


f  (1/s) 


Figure  5.1:  Map  showing  the  experiments  performed  on  the  computer  as  seen  in  N-f 
space;  experiments  in  which  hetons  were  evident  are  denoted  by  the  -f  symbol,  those 
in  which  hetons  were  not  apparent  are  denoted  by  the  o  symbol. 
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